NMR Time Reversal as a Probe of Incipient Turbulent Spin Dynamics 
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We demonstrate time reversal of nuclear spin dynamics in highly magnetized dilute liquid 3 He- 4 He 
mixtures through effective inversion of long-range dipolar interactions. These experiments, which 
involve using magic sandwich NMR pulse sequences to generate spin echoes, probe the spatiotem- 
poral development of turbulent spin dynamics and promise to serve as a versatile tool for the study 
and control of dynamic magnetization instabilities. We also show that a repeated magic sandwich 
pulse sequence can be used to dynamically stabilize modes of nuclear precession that are otherwise 
intrinsically unstable. To date, we have extended the effective precession lifetimes of our magnetized 
samples by more than three orders of magnitude. 
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Elementary treatments of nuclear magnetic resonance (NMR) ignore collective effects. Individual spin dynamics 
are assumed to be independent of the sample magnetization M, and are thus governed by linear differential equations. 
This approximation is justified for many - but certainly not all - practical applications. One of the most well known 
and pervasive counter examples is the phenomenon of radiation damping [l[ , in which the emf induced in the pickup 
coil by M drives a current; the magnetic field associated with this current in turn acts on M, driving it into alignment 
with the static field Bo in a nonlinear manner. A less common but more insidious problem arises when the magnetic 
field produced by the magnetized sample becomes large enough to directly influence spin dynamics @, d, 0, H, [f| . 
This leads to a rich variety of nonlinear effects that range from spectral clustering to precession instabilities and 
spin turbulence. The former arises when small-angle NMR tipping pulses are applied to the sample (generating 
small angular displacements of M from B ), and is manifest by the spontaneous appearance of long-lived geometry- 
dependent modes of coherent nuclear precession. The latter occurs when large-angle tipping pulses are employed, 
and involves an exponential growth in the complexity of spatially inhomogeneous magnetization patterns. Far from 
being mere curiosities, phenomena arising from nonlocal interactions (including the joint action of 'distant dipolar 
fields' and radiation damping @, @|) threaten to limit (or profoundly alter) the operation of state-of-the-art high- 
field high- resolution NMR spectrometers. They are equally important for understanding the dynamics of polarized 
quantum fluids including Bose- Einstein condensates superfluid 3 He EH, degenerate 3 He- 4 He mixtures flij ]. 
and two-dimensional H gases [li|. From yet another perspective, the collective effects induced by distant dipolar 
fields can be used to amplify weak spin precession signals 14] and to enhance the sensitivity of precision searches for 
CP-violating permanent electric dipole moments 0]. 

Here we describe NMR time reversal experiments that - for the first time - shed light on the extent to which the 
deleterious effects of spin turbulence can be controlled or suppressed. Moreover, they lay the foundation for a system- 
atic and quantitative investigation of the onset of dynamical magnetization instabilities under diverse conditions. Our 
experiments are performed on dilute non-degenerate mixtures of hyperpolarized 3 Hc in liquid 4 He. They are based on 
the use of magic sandwich (MS) pulse sequences 15, continuous bursts of tipping pulses that are appropriately 
timed and phased so as to refocus spin dynamics associated with dipolar interactions, leading to the formation of 
echoes. The MS pulse cycle and numerous variants are indispensable tools in the field of solid-state NMR but, to the 
best of our knowledge, have not been used previously for liquid-state NMR 17 1. 

Our samples are confined to a slightly-prolate 0.44 cm 3 spheroidal Pyrex cell treated with Cs metal to suppress 
wall relaxation [l8|, [lj|, and are maintained at temperatures T~ 1 K. The sample cell communicates with a room 
temperature metastability-exchange optical pumping cell and various gas reservoirs via a narrow Pyrex tube. The 
entire apparatus is immersed in a 2.3 mT magnetic field shimmed to ±20 ppm over the sample cell and actively 
stabilized to ±10 ppm against fluctuations in the background laboratory field. We control the volume of liquid 
admitted to the cell (typically 40 % full), the concentration X of 3 He in the mixture (typically 1-5 %), and the 
nuclear polarization (typically < 40 %). We can thus manipulate the shape of the sample, the initial amplitude 
of its magnetization M, and the 3 He diffusion coefficient D. Note that it is convenient to report M in terms of a 
characteristic dipolar field £?dip = HoM or dipolar frequency Fdi p = 7-Bdip/27r, where is the permeability of free 
space and 7 is the 3 He gyromagnetic ratio. 
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Our pulsed NMR experiments are performed at 74 kHz using orthogonal transmit and receive coils. The transmit 
coils are actively-shielded [2(| to suppress the induction of eddy currents in nearby conducting structures. This enables 
us to apply the intense and sustained rf (i.e. Bi) fields that are necessary for NMR time reversal experiments without 
causing the sample to warm. It also ensures that eddy currents do not distort the applied rf field, which is uniform 
to within ±200 ppm over the sample cell. The receive coils are tuned in a tank configuration with a quality factor 
Qo ~ 10. The filling factor r\ is of order 0.02 when the cell is completely full, and an externally adjustable 'Q-spoiling' 
resistor is used to ensure that radiation damping fields are small compared to those generated by the magnetized 
sample. Typically a value Q ~ 2 was employed for the work presented here. 
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FIG. 1: [color online] FID response of polarized 3 He in liquid 4 He (a) to a 90° rf tipping pulse and (b) to a 90°-t-180° sequence. 
The time at which the FID amplitude S drops to half of its initial value So is denoted Ti/ 2 , which is of order 70 ms for these 
examples. Trace (b) has been displaced upward by 0.4 mV for the sake of clarity. The inset shows the deviation of trace (b) 
from its initial amplitude. These data were acquired at T = 1.14 K, with X ~ 2 % and -Bdip ~ 0.9 /xT for trace (a) and 
X ~ 5 % and B dip = 0.7 fiT for trace (b). 



Figure Q] summarizes data from two experiments that illustrate some of the unusual behavior exhibited by highly 
magnetized liquids. The first experiment involves application of a single 90° rf tipping pulse. The corresponding 
free induction decay (FID) is shown as trace (a). The signal amplitude S is roughly constant for several tens of 
ms, during which time the precessing transverse nuclear magnetization is necessarily uniform. This is followed by a 
precipitous drop that culminates in a number of small revivals. Similar experiments performed over a range of initial 
conditions indicate that the characteristic time scale for the abrupt decrease in signal amplitude (i.e. T 1/2 ) is inversely 
proportional to -Bdip 21 1 • 

Evidence that this behavior is not dominated by terms in the equation of motion linear in M is provided by the 
second experiment, in which the initial 90° pulse is followed 70 ms later by a 180° pulse. Comparison of traces (a) and 
(b) reveals that the 180° pulse has little if any effect on the FID envelope. In particular, there is no obvious sign of a 
Hahn (or inhomogeneous) spin echo 140 ms after the start of the experiment. Closer examination of these data shows 
that the deviation of the FID amplitude S from its initial value So grows exponentially with time t (inset to Fig. [T]): 
i.e. 5*o — S ~ exp(2rt). This behavior is expected for an infinite spin system dominated by distant dipolar interactions 
Q , and has been interpreted as the onset of spin turbulence Q . We have performed exact numerical simulations of 
these experiments for small numbers of dipolar-coupled spins diffusing on a cubic lattice (up to 32 x 32 x 32 sites) 
with periodic boundary conditions [211 ] . and observe excellent agreement between measured and calculated growth 
rates. In both cases T/i^ip ~ 2.3, which is somewhat smaller than the value 2y/2n/3 ~ 3.0 predicted for an infinite 
spin system [6|. 
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The experimental challenge presented by these data is to devise a method for probing the spatial and temporal 
evolution of the nuclear magnetization before, during, and after the point at which it spontaneously and abruptly 
averages to zero in Fig. [T] The approach we explore here as a first step toward this goal is the application of MS 
pulse sequences to reverse the effective time evolution of dynamics associated with dipolar interactions, and thereby 
attempt to generate spin echoes. Thorough descriptions of the manner in which time reversal is accomplished can be 
found in the literature [15|. For our purposes the essential features of these multiple-pulse sequences are that (a) the 



applied rf field must be large compared to the dipolar fields and (b) the effective rate at which the dipolar-coupled 
system evolves in time as the rf fields are applied is scaled by a factor of —1/2. 
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FIG. 2: [color online] Spin echo formation: A 90° pulse is applied at t = 0, after which the magnetization is allowed to evolve 
for r = 70 ms. MS sequences of various durations are then used to progressively refocus the magnetization. For traces (a) - (c), 
the amplitude of the FID immediately after the MS probes the rotary echo formed in the doubly-rotating frame. The symbols 
■ , •, a indicate pairings between points on the initial FID to which the state of the magnetization is nominally returned and 
the ensuing free precession signal. For traces (d) - (f), the amplitude of the FID probes the magic sandwich echo formed in 
the rotating frame after the magnetization is first 'unwound' back to t = — r before being 'released' and allowed to undergo 
further free evolution. These data were acquired under conditions similar to those outlined in Fig. [1] Inset: A classic magic 
sandwich [l5l ] consists of a continuous train of (180°,, — 180°,) rotation pairs sandwiched between a leading 90°, and a trailing 
—905) rotation. We have used 360° (rather than 180°) rotations to reduce the accumulation of phase errors associated with 
offsets between the rf and Larmor frequencies. Each 360° rotation is accomplished in 2 ms, corresponding to B\ — 15 /iT. 
Note that for trace (c) the sense of the terminal 90° rotation was reversed (i.e. 90°, instead of —90°,) in order to also refocus 
interactions with a linear dependence on M. This does not significantly alter the initial amplitude of the echo or its subsequent 
growth. 



Figure [5] shows data from several experiments that are analogous to the ineffective 90°-t-180° sequence of Fig. [T] 
except for the fact that the 180° tipping pulse has been replaced with MS sequences of various durations. The rf 
used for these sequences is typically adjusted to be within 1 Hz of the 3 He Larmor frequency. The most striking and 
important feature of these data is that large NMR signals are detected long after the net transverse magnetization 
normally would have averaged to zero. The upper half of Fig. [3 shows data from three experiments that progressively 
'unwind' the state of the dipolar-coupled system back to times t = 2r/3, t — t/3, and t — 0, as indicated by the 
paired-symbols ■, •, and a. The initial amplitudes of these echoes nominally probe the amplitude of the 'rotary echo' 
[l5l [22} that is formed in the doubly-rotating frame [23| during the MS sequence. In the first case - trace (a) - the 
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magnitude of the transverse magnetization is returned to a value that is close to that which is anticipated. Moreover, 
the shape of the ensuing FID mimics the original decay (cf. Fig. []}. In the last case - trace (c) - the amplitude of 
the signal following the MS starts out noticeably smaller than anticipated but then unexpectedly grows rather than 
decaying. Clearly the deficit in the initial echo amplitude is not associated with an irreversible process. The shape of 
the FID labeled (b) is qualitatively intermediate between traces (a) and (c). The lower part of Fig. [5] shows data from 
three experiments that nominally 'unwind' the state of the dipolar-coupled system back to a time t = — r, resulting 
in the formation of a 'magic sandwich echo' in the rotating frame (TBj after a further free evolution period of duration 
t. Increasing the initial magnetization increases the rate at which magnetization patterns develop, and so this series 
of experiments effectively probes the sample at successively later stages of evolution. 

It is not surprising that we observe imperfect refocussing of the transverse nuclear magnetization in these exper- 
iments. Irreversible processes (e.g. diffusion) necessarily compete with the spatiotemporal evolution of the magne- 
tization and progressively attenuate structures on larger and larger scales as time progresses. At the same time, 
while the MS tends to refocus the family of magnetization patterns that develop during the initial FID, one expects 
a complementary family of instabilities [6( to be established and grow during the time reversal sequence. Technical 
imperfections in the time reversal sequence will also contribute to echo attenuation. A detailed and systematic study 
of these rich and complex effects is required before quantitative statements can be made regarding the onset and 
development of spin turbulence. 

The fact that spin echoes can be generated in highly magnetized liquids leads one to contemplate the use of more 
sophisticated NMR tools. To illustrate this point, we describe here a simple but significant extension of the basic 
MS sequence that clearly merits further investigation. A conventional CPMG sequence (90 o -t-180°-2t-180°-2t . . .) 
acts on a system dominated by interactions that are linear in M, repeatedly inverting the phase of the precessmg 
magnetization to produce an echo train. By analogy, a Repeated Magic Sandwich (RMS) sequence consisting of a 
train of MS cycles interspersed with appropriate delays should likewise act on a system dominated by distant dipolar 
interactions to form an echo train. The validity of this conjecture is borne out by the data shown in Fig.[3l The RMS 
sequence used to generate these data is of the form (90°-t-MS-2t-MS-2t . . .), where the duration of each MS cycle is 
4r (i.e. long enough to effectively 'reverse' the time evolution of the dipolar coupled system for a time 2r). The four 
traces correspond to experiments performed with different delay times r. In each case, the RMS sequence is applied 
for ~ 9 s, at which point M is allowed to evolve without further interruption. 
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FIG. 3: [color online] Spin echo trains: Traces (a) - (d) correspond to echo trains that are observed when RMS sequences with 
delay times r = 16, 24, 32, and 40 ms, respectively, are applied to the magnetized liquid. The MS cycles used to obtain these 
data are identical to those described in the caption of Fig. [2] except for the facts that each one terminates with an additional 
1805) rotation (cf. trace (c) of Fig. to refocus interactions linear in M, and that the symmetry of alternate MS cycles is 
inverted to reduce the accumulation of phase errors. That is, odd MS cycles are of the form [90)), n x (360J, —360^), 90$)] while 
even MS cycles are of the form [—905), n x (—360)), 360)))), —90))]. These data were acquired at T ~ 1.1 K, with X ~ 4.5 % and 
-Bdip ~ 0.8 fiT. In the limit B d i P -> 0, CPMG-based measurements yield D~2x 10" 7 m 2 /s. 
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Figure [3] reveals two distinct phases in the evolution of the magnetization during the application of RMS sequences. 
Each echo train starts with a period that we associate with imperfect refocussing of the magnetization. That is, 
behavior similar to the examples presented in Fig. O Increasing the length of time during which the magnetization 
is allowed to evolve between successive magic sandwiches results in greater losses. This initial period is followed by a 
second phase during which magnetization losses continue to be incurred, but at significantly lower rates. Comparison 
with Figs. [1] and [5] reveals that the apparent lifetimes of these dynamically sustained precession signals are in some 
cases more than 3 orders of magnitude longer than the characteristic coherence time (e.g. T 1/2 ) of the freely evolving 
magnetization. Exact numerical simulations of these experiments (analogous to those described above) have been 
used to investigate the influence of diffusion and rf field amplitude. They reveal similar - but not identical - behavior. 
The initial period during which losses are incurred is observed, as is the transition to a second evolutionary phase 
and the dependence of this transition on both r and -Bdip- On the other hand, our model calculations yield a 
steady state solution for which we are not able to resolve a decay rate: that is, a plateau. This suggests that our 
experimentally-determined echo trains may yet be limited by systematic rather than intrinsic effects. One possibility 
is the small but finite heat load placed on our cryostat by eddy currents induced during RMS sequences, which would 
in turn lead to evaporation of 3 He atoms. Another possibility is imperfect refocussing of the magnetization caused 
by imperfect phasing of rf tipping pulses. Finally, it is worth commenting on the terminal FID at the end of each 
RMS sequence. Analysis of these decays yields values of T that are proportional to the terminal echo amplitude 
and consistent with those observed following a single 90° tipping pulse (cf. Fig. Q}. This strongly suggests that the 
precessing magnetization is highly uniform at the peak of each revival, and that the 'seeds' from which instabilities 
ultimately develop @ are small. 

In summary, we have demonstrated the production of NMR spin echoes and spin echo trains in magnetized liquids 
under conditions where nuclear spin dynamics (in the rotating frame) are dominated by nonlinear and nonlocal 
(or distant) dipolar interactions. Our experiments, which are based on the use of magic sandwich time reversal 
pulse sequences, probe the rich spatiotemporal evolution of the nuclear magnetization distribution that results from 
competition between nonlinear spin dynamics and diffusion. They promise to yield quantitative information about 
the progression and extent of disorder, as well as new methods to control or suppress the onset of instabilities. 
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